Behaviour of muonium in synthetic diamond 
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1 Introduction 

Diamond with its unsurpassed mechanical strength, thermal conductivity, and radia- 
tion hardness is a promising material for radiation detectors, for electronic and optical- 
electronic components able to withstand high heat and radiation loads. Great advances 
have been made over the last years in technology of manufacturing synthetic single-crystal 
diamond and diamond films [HE]. The properties of synthetic diamond should be com- 
qq ■ prehensively studied to find out the scope of its practical application. 

Man-made bulk diamond samples (single crystals and films) contain about (or more 
than) 10 17 cm -3 dissolved hydrogen atoms. Recently, surface conductivity was observed 
in a diamond sample saturated in a hydrogen atmosphere j3]. However, conventional 
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methods do not yield sufficient information on the behaviour of hydrogen in diamond. 

Polarized positive muons jjl + can be used to imitate and study the behaviour of an 
hydrogen atom in matter (see, for example, [HE]). A positive muon in matter may 
pick up an electron and form muonium (Mu = /i + e~). Since the muon mass is about 
1/9 than that of the proton, muonium can be considered as a light isotope of hydrogen. 
Theoretical calculations [HIE] show that in semiconductors with diamond crystal structure 
the muonium may occupy tetrahedral and octahedral interstitial sites or is localized in 
the middle of the axis between two host atoms (bond-centre site). It was experimentally 
found (see [8j[9]) that in diamond the muonium occupies the tetrahedral site (commonly 
referred to as normal muonium and designated as Mu or Mur) and the bond-centre site 
(commonly referred to as anomalous muonium and designated as Mu* or Mubc)- One 
more /iSR signal corresponding to the diamagnetic state of /i + was observed in experiments 
(yU + -state). 

In diamond Mut possesses isotropic symmetry with a hyperfine constant A-^/h = 
3711 ± 21 MHz [8] (A hi /h = 3693 ± 83 MHz [ID]), which is slightly smaller than that 



for muonium in vacuum Aq/H = 4463302.765 ± 0.053 kHz [TTJ. At the bond-centre site 
the hyperfine interaction is anisotropic with axial symmetry relative to the [111] axes of 
the crystal. The hyperfine interaction constants of Mubc are A\Jh = 167.98 ± 0.06 MHz 
and A±/h = —392.59 ± 0.06 MHz [§]. For the muonium bond-centre site the muon spin 
precession frequency in a magnetic field depends on the angle between the magnetic field 
and the [111] axes of the crystal. However, theoretical calculations [6j[7] reveal that for 
Mubc there is a certain 'magic' field where the muon spin precession frequency is nearly 
independent on the orientation of the crystal. This effect can be used to observe a Mubc 
signal at the 'magic' field in a powder sample. 

The above-mentioned values of the hyperfine constants of the muonium were obtained 
for two natural diamond samples: la-type single-crystal diamond and Ila-type diamond 
powder with grain size 1-6 /zm. In the present work the behaviour of the positive muon 
in two synthetic diamonds is studied. 



2 Measurements 

The measurements were carried out at the GPS and Dolly spectrometers located in the 
7rM3.2 and 7rEl muon beams of the Paul Scherrer Institut (PSI, Switzerland). In the 
7rM3.2 beam line the muons were transversely polarized — the angle between the muon 
spin and the muon momentum was ~ 70°. 
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Figure 1: Location of the apparatus in the PSI muon beams. Up, Dw, Fw and Bw 
are the up, down, forward and backward telescopes for detection of positrons. M is the 
counter for registration of the incoming muons, H is the magnetic field, S is the sample, 
Sfj_ indicates the muon spin polarization direction, a) At the GPS spectrometer the muon 
beam is transversely polarized and the magnetic field is parallel to the incoming muon 
beam momentum, b) At the Dolly spectrometer the muon beam longitudinally polarized 
and the magnetic field is perpendicular to the picture plane. 



The location of the apparatus in the PSI muon beams is schematically shown in FigJTJ 
Muons with energy ~ 4 MeV pass through the hole in the Bw counters and a thin M- 
counter and stop in the sample S. The positrons from the muon decay were detected by 
the Bw (backward), Fw (forward), Up(up) and Dw(down) telescopes located around the 
sample (see Fig. [I]). A magnetic field transverse to the muon spin polarization at the 
sample was created by Helmholtz coils. The long-time stability of the magnetic field was 
better than 10 -4 . 

Measurements were carried out using a helium flow cryostat that allowed changing 
temperature of the sample within the range of 4.2-300 K. The temperature of the sample 
was stabilized with an accuracy better than 0.1 K. Measurements were carried out in 
'magic' magnetic field Hm = 0.4073 T at the GPS spectrometer and in H = 7.5 mT at 
the Dolly spectrometer. 

In these measurements two samples (designated as D3 and D5) were used. Sample 
D3 with the total mass of 2 g was composed of a few Ila-type single-crystal diamonds. 
Orientation of the crystals in the sample was random (MSC - mosaic single crystal). 
The concentration of nitrogen in sample D3 was about 0.5 ppm (in diamond, 1 ppm 
corresponds to 1.76 • 10 17 cm -3 ). The crystals contain few wt% of Ni as isolated inclusions 
visible by eye. The diamond crystals with dimension 1-4 mm were synthesized at the 
Technological Institute for Superhard and Novel Carbon Materials (Troitsk, Russia) at 
high pressure and high temperature with synthetic diamond as seed material [12J. Sample 
D5 (CVD) of total mass 2 g was a batch of pieces of several diamond films 0.5 mm thick 
and about 10 x 5 mm in lateral size. The diamond films were produced by a microwave 
plasma assisted CVD technique in CH 4 /H 2 mixtures [T3]. The main impurities in sample 
D5 were hydrogen (65 ppm) and nitrogen (1.5 ppm). Other impurities were less than 
0.1 ppm. The crystalline axis [110] perpendicular to the film plane was the predominant 
grain orientation. 

The evolution of the muon polarization P(t) was studied by measuring the time dis- 
tribution of positrons from the decay of muons /i + — > e + + v e + v^ . In the general case the 
time distribution of the positrons (with respect to muon stop in sample) can be presented 
by the function 

N(t) = N ex P (-t/r^)(l + a/3 -G(t)) + Bg, (1) 

where N is proportional to the number of muons stopped in the sample, r M is the muon 
lifetime, G(t) is the polarization of the muon at the decay, a depends on the parameters 
of the /iSR setup and is close to 1, Bg is background. 

The function G(t) depends on the experimental condition: for measurements at H = 
7.5 mT and H = 0.4073 T the respective explicit expressions for G(t) are 

G(t) = P„(0) exp(-iy) cos(27rz/ M t + <p ) 

+ P T (0)exp(-R T t)cos(2nQ 1 t + ip 1 )cos(2nQ 2 t + ^2), 

G(t) = P M (0) exp(-iy) cos(2vrz/ M t + ip ) + P BC (0) exp(-R BC t) cos(2iTu M t + </?i), (3) 

where Qi = (^23 + ^i2)/2; ^2 = (^23 — ^i2)/2; v ik is the frequency of the transition between 
levels of hyperfine structure of Mur; variables with indices /x, T and BC refer to the muon 



in the diamagnetic, normal muonium and anomalous muonium states respectively; -P?,(0) 
is the muon polarization at t — 0; Ri is the relaxation rate of the muon spin; vyi is the 
muon spin precession frequency for the Mubc fraction in the 'magic' magnetic field; (pi is 
the initial phase of the muon spin precession. 

3 Results and discussion 

Figure [2] shows the results of Fourier transformation of the /iSR spectra for sample D5 
measured in the magnetic field 7.5 mT at the temperatures 50 and 250 K. The peaks 
at ~ 100.4 and ~ 105.8 MHz correspond to the transitions between the levels of the 
hyperfine structure of the muonium Mur- The narrowing of the lines z/12 and ^23 with 
increasing temperature was observed for both samples D3 and D5. 
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Figure 2: The result of the Fourier analysis of the /xSR histograms for diamond film D5 
measured in 7.5 mT magnetic field at the temperatures a) 50 and b) 250 K. 



sample 


50 K 100 K 200 K 


250 K 


D3 (MSC) 


19.1 ±0.8 


13.4 ±0.9 


D5 (CVD) 


12.4 ±0.6 17.3 ±1.4 12.3 ±0.6 


7.7 ±0.6 



Table 1: Muon spin relaxation rate Rt (MHz) for the Mut fraction in synthetic diamond 
at different temperatures. 



The /iSR spectra measured in the magnetic field 7.5 mT were fitted by (CQ) with 
polarization function (J2J) and values of the parameters -P^(O), -Pt(O), -Rt, ^, ^i and f2 2 
were found (see Table [1] and Table [2]). It was found that R^ is equal to zero, and final 
values for other parameters were obtained by fitting the experimental data with fixed 
Ru = 0. The values f2i and f2 2 presented in Table [2] correspond to more precise results 



obtained at 250 K. Hyperfine constants were calculated as 
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sample 


v» MHz 


n u MHz 


fi 2 , MHz 


A M , MHz 


D3 (MSC) 


0.996 ±0.005 


103.73 ±0.12 


2.95 ±0.22 


3715 ± 277 


D5 (CVD) 


0.992 ±0.004 


103.13 ±0.06 


2.75 ±0.10 


3940 ± 143 



Table 2: Muon spin precession frequency and hyperfine interaction constant A^ for Mux- 

The values of the hyperfine constant for Mut in samples D3 and D5 are smaller than 
for muonium in vacuum, and they are in agreement with the value A^ = 3711 ± 21 MHz 
for natural diamond powder [8] within the accuracy of the measurements. However, a 
significant difference was observed in the behaviour of muonium in synthetic diamond in 
comparison with Ila-type natural powdered diamond [8]: 1) the muon spin relaxation rate 
i?T for samples D3 and D5 is approximately ten times higher than for the natural sample; 
2) in distinction to results [8], Mut was observed at temperatures higher than 150 K, and 
at 250 K the relaxation rate decreases slightly (see also Fig. [2]). At the same time the 
present results for the relaxation rate in synthetic samples D3 and D5 are close to those 
observed in p3] for Ha- and lib-type natural diamonds. 
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Figure 3: Muon spin precession frequencies observed in synthetic diamond D3 in the 
'magic' magnetic field 0.4073 T at 100 K. 



Figure |3] shows the results of the Fourier analysis of the /iSR-data measured by the 
Up telescope of the GPS spectrometer for mosaic single crystal diamond sample D3 in 
the 'magic' magnetic field 0.4073 T at 100 K. There are two peaks: one corresponds to 
the precession of the muon spin in the diamagnetic state and the other is due to the 



precession of the muon spin in the anomalous muonium in the 'magic' magnetic field. 
Muon spin precession frequencies of Mut too high to be observed in the "magic" field. 
The /iSR histograms collected by the Up and Dw telescopes were fitted with polarization 
function ([3]) and values of the parameters -Rbc> ^ anc ^ ^m were extracted: Rbc = 2.1 ± 
0.5 MHz, V/l = 55.224 ± 0.005 MHz and v u = 142.8 ± 0.1 MHz. As for the normal 
muonium in synthetic diamond D3, the relaxation rate of the muon spin in anomalous 
muonium is an order of magnitude higher than that observed earlier in Ila-type natural 
powder diamond (0.25 ± 0.10 [8]). 




Figure 4: The time dependence of the muon polarization in diamond D3 in the magnetic 
field 0.4073 T for the case of longitudinal field measurements (LF). The data are corrected 
for the muon life time r^ = 2197 ns. 

Since in the 7rM3.2 beam line the angle between the muon spin and the muon mo- 
mentum is ~ 70° and the external magnetic field is parallel to the muon momentum, the 
yuSR-histograms collected by the Fw and Bw telescopes (see Fig. [T^, ) contain information 
about evolution of the longitudinal (with respect to the magnetic field) component of 
muon polarization (LF measurement). In Fig. H]the /iSR histogram for sample D3 (MSC) 
collectedted by the Fw telescope at the magnetic field H = 0.4073 T and T = 100 K are 
presented on two different time scales. There are two components of muon polarization 
in the histogram: oscillating and non-oscillating. The Fourier analysis reveals only one 
frequency in the Fw and Bw /iSR histograms which corresponds to muon spin precession 
of Mubc i n the 'magic' magnetic field. It is evident that the oscillating component of 
polarization is due to the Mubc fraction. At the same time both Mut and Mubc frac- 
tions may contribute to the non-oscillating component of the Fw and Bw /iSR histograms. 
These histograms were fitted by the polarization function 

G(t) = P(0) expH/Tx) + Pbc(0) exp(-t/T 1 BC )) cos(2vri/ M t + <f), (5) 

where 7\ and T® c are the longitudinal relaxation times for the non-oscillating and oscil- 
lating components of muon spin polarization. 

The following values of the muon spin polarization parameters were found: 1/Ti = 
1.07 ± 0.02 MHz, l/7f c = 1.14 ± 0.14 MHz, u M = 142.68 ± 0.06 MHz. As was expected, 
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Table 3: Fractions of muon states (/i + , Mux, Mubc) an d missing fraction (MF) observed 
in samples of natural (n) and synthetic (s) diamonds. SC is single crystal diamond, and 
MSC is a sample composed from few single-crystal, CVD is chemical vapor deposited 
diamond film. *-we estimated according to the data presented in original publication. 
**-present data. 



the value of the z/m obtained from LF measurements is in good agreement with the results 
of the TF measurements (data collected by the Up and Dw telescopes). 

To find the fraction of muon formation in the /i, Mux and Mubc states in diamond 
the values obtained for a ■ -P M (0), a ■ -Pt(O) and a ■ -Pbc(O) were normalized to the muon 
polarization for the silver sample measured under the corresponding experimental con- 
dition. It was found that for the Up, Dw and Fw telescopes of the Dolly spectrometer 
a ■ PAg(0)/3 = 0.270 ± 0.001 and for the Up and Dw telescopes of the GPS spectrometer 
a ■ PAg(0)/3 = 0.240 ± 0.001. In the calculation it was taken into account that in the 
transverse magnetic field only 0.5 of the muon polarization in Mur and in 'magic' mag- 
netic field 11/30 of the muon polarization in Mubc [8] gives experimentally observable 
precession signal. The fractions of muons formed in the fi + , Mur and Mubc states in 
diamond samples D3 and D5 are presented in Table [3] in comparison with the analogous 
data for various diamond samples studied earlier. 

The data in Table |3] are divided into the three sets. The first set includes data for the 
la-type single-crystal and Ila-type powdered natural diamonds, the second set presents 
data for the Ila-type and lib-type single-crystal natural diamonds and the third set is for 
the synthetic diamond. The Mux fraction is minimal and the missing fraction is maximal 
for the first set of samples. For the second set of samples the missing fraction of muon 
polarization is close to zero. In the synthetic samples the Mubc fraction is 2-3 times 
smaller than in the Ha- and lib-type single-crystal natural diamonds. The large missing 
fraction in the Ila-type powder sample (with grain size 1-6 /zm) was explained by diffusion 
of the Mut to the crystalline boundaries [7]. 



It was experimentally shown that a large missing fraction of the muon polarization 
in the la-type diamond is due to one more muonium state in the nitrogen-rich dia- 
mond [HUH?]. As is well known, la-type diamonds contain about 0.1% (1000 ppm) 
nitrogen impurities and they are aggregated to large defects forming A- and B-centres. 
At A-centres two nitrogen atoms occupied two nearest substitution sites. The calcula- 
tion given in [IT] shows that the muonium may occupy the bond centre site between the 
neighboring nitrogen atoms in the case of A-type defects and the site between the nitrogen 
molecules in the case of the B-type defects. 

The concentration of nitrogen atoms in the Ha- and lib-type natural diamond is 
approximately three orders of magnitude lower than in la-type diamond. The fact that 
the missing fraction of muon polarization in the Ha- and lib-type natural diamond is close 
to zero indicates that nitrogen does not play any noticeable role in formation of muonium 
fractions in these samples. 

In the synthetic diamond the concentration of nitrogen is close to that in the Ilb- 
and Ila-type natural diamond. Therefore it is reasonable to compare the results for the 
samples studied in this work with those for Ha- and lib-type natural diamond. As is 
seen from Table |3j in the synthetic samples the Mup fraction is comparable, the Mubc 
fraction is 2-3 time smaller and the missing fraction is a few times larger than in the Ila- 
and lib-type single-crystal diamonds. The fact that the missing fraction in the synthetic 
diamond is 2-3 times smaller than in the nitrogen-reach la-type natural sample and 2-3 
times larger than in the Ha- and lib-type samples indicates presence of some other defects 
in the synthetic diamond in comparison with the natural one. The observed difference 
may be due to the high concentration (more than 10 ppm) of incorporated hydrogen and 
due to a small crystalline size (a large effective surface) in the CVD films and due to the 
presence of Ni includes in sample D3 (MSC). 

4 Summary 

It is found that the hyperfine constant for interaction of the muon and the electron 
magnetic moments on Mup and Mubc i n the synthetic diamond is close to those in the 
natural diamond. The longitudinal relaxation rate of the muon spin in the Mubc state is 
measured for the first time: 1/Tf c = 1.07 ± 0.02 MHz. The Mut fraction of the muon 
polarization and the relaxation rate R<r for the synthetic samples are comparable with 
those for the Ha- and lib-type natural diamond. 

The fact that 1) the missing fraction of the muon polarization in the synthetic samples 
is 2-3 times smaller than in the la-type natural sample and 2) the Mut fraction and the 
relaxation rate Rt for the synthetic samples is comparable with those for the Ha- and 
lib-type natural diamond indicates to the low concentration of such defects as A- and 
B-centers in synthetic samples. 

We are grateful to the Directorate of the Paul Scherrer Institute for allowing us to 
carry out the experiments at PSI. 
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